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The MiniBooNE collaboration has recently measured the neutrino charged current (CC) 
quasielastic (QE) cross section on ^^C for a neutrino beam with average energy of 788 MeV [1]. 
The comparison of these results with a prediction based on the relativistic Fermi gas model using 
the standard value of the axial cut-off mass Ma = 1-03 GeV/c^ reveals a substantial discrepancy. 
In the same model the larger value of Ma = 1-35 GeV/c^ is needed to account for data. For 
nuclear physicists accustomed to the complexity of the nuclear system, this anomaly is likely to 
reflect the many-body aspect of the problem. Indeed Martini et al. have pointed out [2j that, 
depending on the detection method, certain types of inelastic events can simulate quasielastic 
ones. As the MiniBooNE experimental set-up defines a "quasielastic" event as one in which 
only a muon is detected, the ejection of a single nucleon (a genuine quasielastic event) is only 
one possibility, and one must in addition consider events involving a correlated nucleon pair 
from which the partner nucleon is also ejected. This leads to the excitation of 2 particle-2 hole 
(2p-2h) states; 3p-3h excitations are also possible. The inclusion in the quasielastic cross section 
of events in which several nucleons are ejected, leads to an excess over the genuine quasielastic 
value. Martini et al. have argued that this is the likely explanation of the anomaly showing 
that their evaluation can account for the excess in the cross section without any modification 
of the axial mass. After this suggestion, a number of articles have discussed the problem of 
multinucleon emission and whether it could account for the anomaly. The aim of this work is 
to review the actual theoretical status on this subject. 

The theoretical studies of 2p-2h in connection with CCQE at MiniBooNE kinematics are 
actually performed essentially by three groups. There are the works of Martini et al. [21 [3l H], 
the ones of Amaro et al. O [6| and the ones of Nieves et al. [3 |8]. Recently it appears also 
a paper of Bodek et al. [9] considering a more phenomenological approach related to electron 
scattering. Two-nucleon emission via two-body currents is not a new idea related to neutrino- 
nucleus interaction, it was studied since the eighties in connection to the electron scattering, 
photon and pion absorption. It is for example well known that the absorption of pions by nuclei 
is a two-nucleon mechanism. Concerning electron scattering, several theoretical calculations 
have shown that 2p-2h are needed to reproduce the inclusive cross section. In particular they 



contribute to fill the so called "dip" region between the quasielastic and the delta peaks. They 
affect the isospin spin-transverse (i.e. the magnetic) response. An exhaustive bibliography 
on these studies is beyond the purposes of this work. In the following we just mention the 
sources and the references of the calculations of the three groups mentioned above. The np- 
nh contributions in the papers of Martini et al. are obtained starting from the microscopic 
calculations of pion absorption at threshold [lOJ, of pion and photon absorption [llj and of 
the transverse response in electron scattering [12j. They give two different parametrization of 
this np-nh contribution (which leads to practically identical neutrino total "quasielastic" cross 
section), one more relates to [lUj . the other to [T5]. In the following we present the results 
just obtained with this second parametrization which is characterized by a more satisfactory 
energy and momentum transfer dependence. We refer to ^ for details. The 2p-2h contributions 
considered by Amaro et al. are taken from the full relativistic model of [13] related to the 
electromagnetic responses. The work of Nieves et al. can be considered as a generalization of 
[14], developed for the electron scattering, to the neutrino scattering. The contributions related 
to the A-MEC excitations are taken, as in the case of Martini et al. from [11]. Considering 
these three different models, it exists some differences already of level of genuine QE which 
can be particularly important when one compares the double differential cross sections. Amaro 
et al. considered the relativistic superscaling approach based on the superscaling behavior 
exhibited by electron scattering data. The models of Martini et al. and Nieves et al. are 
more similar: they start from a local Fermi gas picture of the nucleus. They consider medium 
polarization and collective effects through the random phase approximation (RPA) including 
A-hole degrees of freedom, vr and p meson exchange and g' Landau-Migdal parameters in the 
effective p — h interaction. The treatment of Nieves et al. is relativistic while the original 
one of Martini et al. was non relativistic. Recently it has been improved through relativistic 
corrections Concerning the 2p-2h sector we remind that there exist several contributions 
to two-body currents, see for example [3 [12l [13]. There are the so called pion-in-flight term, 
the contact term and the A-intermediate state or A-MEC term. At level of terminology, some 
authors refer just to the first two terms as Meson Exchange Currents contributions, like in 
|2], but actually the most current convention consists into include the A-term into MEC. In 
order to preserve the gauge invariance of the theory also the nucleon-nucleon (NN) correlation 
contributions must be taken into account. In an infinite system they give rise to divergences. 
The different prescriptions to regularize them lead to a model dependence of the results. In this 
connection Amaro et al. consider only the MEC contributions and not the correlations and the 
correlations-MEC interference terms. Correlations and interference are present both in Martini 
et al. and Nieves et al. even if Martini et al. consider only the A-MEC. On the other hand the 
treatment of Amaro et al. is fully relativistic as well as the one of Nieves et al. while the results 
of Martini et al. are related to a non-relativistic reduction of the two-body currents. Amaro et 
al. consider the contribution only in the vector sector while Martini et al. and Nieves et al. also 
in the axial one. The virtual meson exchanged between the two nucleons is just the pion for 
Amaro et al. while Martini et al. and Nieves et al. includes also g' Landau-Migdal parameters. 
Nieves et al. consider also the p meson exchange contribution. 

Before to review the results obtained in neutrino scattering, let's analyze the transverse 
response in electron scattering. In Fig[T] we present the transverse (e, e') response on carbon 
calculated in the model of Martini et al. for some values of the momentum transfer q as a 
function of the energy transfer uj together with the experimental results \15\ 116] . Note that 
there are already some differences between the two data sets. Concerning the results at g = 300 
MeV/c one can compare the Martini et al. 2p-2h contribution with the one of Fig. 41 of |14]. 
One can observe that they are of the same order of magnitude. At higher q np-nh seems to be 
needed in order to reproduce data in particular in the dip region. Bodek et al. p5] parametrized 
this enhancement in terms of a correction to the magnetic form factor. 
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Figure 1. (e,e') transverse responses of ^^C. Experimental data are taken from [15j and [16]. 
The theoretical calculations are performed in the model of Martini et al. 

Turning to neutrino scattering, FiglJ] shows the "quasielastic" and the inclusive CC cross 
section as measured by MiniBooNE ^ and SciBooNE [I7| respectively. The theoretical results 
of Martini et al. are also reported. As already mentioned the inclusion of the multinucleon 
component leads to an agreement with the "quasielastic" MiniBooNE data. Also the inclusive 
cross section, which is less sensitive to the background subtraction with respect to exclusive 
channels, can be reproduced, at least up to E'jy ~ 1 GeV once the np-nh channel is added. 
Similar conclusions both for QE and inclusive cross section have been obtained by Nieves et 
al. [7]. Concerning the MiniBooNE "QE" total cross section, it can actually be reproduced 
without any change in the axial mass also by Bodek et al. which on the other hand agrees also 
with NOMAD data. Also for Amaro et al. the agreement with the total "QE" cross section 
improves, at least at high neutrino energy, adding to the genuine QE in the SuSA approach 
the MEC contribution. Nevertheless the shape of this cross section is somewhat different with 
respect to the other calculations, maybe owing to the absence of NN correlations contributions. 
The total cross section has been analyzed also in the antineutrino mode by Martini et a/.[3j, 
Nieves et aL[7] and Bodek et a/. [9]. The effect of np-nh as compared to one nucleon ejection 
should be somewhat less important for antineutrino because the spin isospin response which is 
affected by the 2p-2h piece has less weight for antineutrino owing to the negative axial-vector 
interference. This effect is evident in Martini et al. and Bodek et al, less in Nieves et al. Let's 
turn to the MiniBooNE flux averaged "QE" double differential cross section. This is a directly 
measured quantity, hence free from the uncertainty of neutrino energy reconstruction which is 
a model dependent procedure. Theoretical calculations on 2p-2h contributions to this cross 
section have been performed by Amaro et a/.[5l|6], Nieves et al. [8] and Martini et al. [4J. Two 
examples of results of Martini et al. are shown in Fig. [3j For a complete comparison in all the 
cos6 and range measured by MiniBooNE see [1]. The agreement with data is good in all 
this range once the multinucleon component is incorporated. It confirm that it is possible to 
account for the MiniBooNE results without any modification of the axial mass. Similar results 
have already been obtained by Nieves et al. [8\. Although slightly lower with respect to the 
Martini et al. calculations they are also fully compatible with MiniBooNE using the standard 
value of the axial mass. In some kinematical regions the agreement with MiniBooNE arises 
in Martini et al, as well as in Nieves et al., from the opposite actions of the RPA quenching 
and the enhancement from the np-nh contribution. For Amaro et a/.[5il6| the inclusion of 2p- 
2h MEC tends to improve the agreement with data at low angles, but it is not sufficient to 
account for discrepancy at higher angles. The absence of NN correlation contributions might 





Figure 2. Inclusive and "quasielastic" cross 
sections measured by SciBooNE [T7] and 
MiniBooNE [1] compared to Martini et al. 
calculations [2]. 



Figure 3. MiniBooNE [Ij flux-averaged 
double differential QE cross section for 
some values of and cos^. Theoretical 
calculations are the ones of Martini et al. [3] . 



be responsible for residual disagreement. A treatment of these contributions performed by the 
same group in the electron scattering sector [18] seems to suggest that they are of the same order 
of magnitude as the one of 2p-2h MEC and that they are supposed to have a greater impact 
for large neutrino- muon angles. A further enhancement of the cross section expected from these 
contributions would lead to an agreement with MiniBooNE. 

In conclusion all the actual theoretical models treating the quasielastic with the inclusion of 
the multinucleon emission channel suggest that the MiniBooNE results can be explained without 
any modification of the axial mass. The complexity of the nuclear dynamics needed to reproduce 
the data should be taken into account in the algorithm to reconstruct the neutrino energy which 
is crucial in neutrino oscillation experiments. 
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